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re-analyses is discussed, focussing on the complementarity of the two approaches. 
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1. Introduction and Background 

The rif = 5 QCD coupling in the MS scheme at the conventionally defined reference scale 
H = Mz represents one of the fundamental parameters of the Standard Model. The central value 
of the 2008 PDG assessment a.v(Mz) = 0.1176(20), remains strongly influenced by the high- 
precision lattice result, as{Mz) = 0.1170(12), obtained in Ref. [||] from an analysis of UV-sensitive 
lattice observables using the MILC a ~ 0.18, 0.12, and 0.09 fm ensembles. In the last year, a 
number of independent determinations have appeared, from a number of different sources, yielding 
typically somewhat higher values. Specifically, 

• updates of the global EW fit, taking into account the new 5-loop result for the dimension 
OPE contributions yield a,(Mz) = 0.1191(27) [|,|]; 

• the most recent hadronic T decay determinations |^, |7|], whose ranges are encompassed by 
results of Ref. yield a, (Mz) = 0.1 187(16) (~ 2a lower than the earlier result of Ref. ||], 
for the reasons discussed in Refs. ^); 

• the determination based on 2 to 10.6 GeV Rhad values yields as{Mz) = 0.1190 (Ino) 

• an update of the determination from rffi"])!!^^ yields a, (Mz) = 0. 1 190 ( tso) [0 ; 

• recent determinations, a,(Mz) = 0.1198(32) [|lO|], 0.1182(45) 0.1240(33) 10] and 
0.1172(22) [pj|], using shape observables in DIS and e^e^ hadrons, if averaged naively, 
yield a combined result av(Mz) = 0.1193(15). 

In view of these results, and, in addition, the availability of the new MILC ensembles with a ~ 0.15 
and 0.06 fm, it is timely to revisit the earlier lattice analysis. We do so by focussing on three 
observables, log(Wii), log(Wi2) and log(Wi2/Mo) which can be convincingly argued to receive 
only small non-perturbative contributions at the scales of the lattices employed in our analysis. In 
what follows, we first briefly outline the basics of the method employed in Ref. [^], then describe 
our implementation of this approach, and finally present our results. We also discuss briefly the 
differences (and complementarity) between our implementation and that of the other recent similar. 



but not identical, reanalysis by HPQCD [|14|]. An expanded discussion of the work reported here 
may be found in Ref. [|T5|]. 

The authors of Ref. [^] extracted as{Mz) by studying a large number of UV-sensitive lattice 
observables, including the three, log(Wii), log(Wi2) and log(Wi2/Mo) o^i which we focus below. 
The perturbative expansion for such an observable. Ok, is written in the form 

= £ ^/«a(GO'' = £»,«(&) £ ci^am" d-D 

ik) 

with Cq =\, Qk = dk/a the Brodsky-Lepage-Mackenzie (BLM) scale for the observable O^, and 
a any coupling having the same expansion to 0{a^) (with as the usual MS coupling) as the usual 
heavy quark potential coupling, denoted ay below. The coefficients c?i 23 (equivalently, DkiC\ , 



and C2 ) have been computed in 3-loop lattice perturbation theory [ 16 1 for a number of such observ 



ables and, with the corresponding du, tabulated in Refs. |Q, 14, T^]. They are common for all such 
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couplings a. The couplings a also share common values for the first three j5 function coefficients, 
j3o = 9/4, j3i = 4 and p2 = 33.969, where in our normahzation jx^ da{ix) / djl^ = -T.n=oPna"{lJ.) 
with d = a/n. When the expansion of a is further specified to 0{af), jSs is also determined from 
the known values of the 4-loop MS j8 function coefficients, j3o, • • ■ , jSs JXt]]. The a{Qk) appearing 
in Eq. ( |1.1[ ) are determined by the value at a single reference scale, the reference scale value serv- 
ing as a fit parameter for the analysis. It is, of course, important to remove any non-perturbative 



contributions to the observable in question in order to make use of Eq. (|L1|). The reliability of the 
analysis will be greatest when such non-perturbative subtractions are small. 

Somewhat different choices for a are made in the two recent reanalyses. We denote our choice 
by aT, and that of Ref. |[T^] by ay. The relation between tty and as, to 0{a^), is of the form [|18|] 

aP{q^) = a^ipL") [1 + JCi(AtV9')«s(M') + K:2(AtV'?')«.(At')] (1-2) 



with the expressions for Ki^ix) given in Ref. |[T8|]. The tif = 3 version of the RHS of Eq. (1.2), 
with jj.^ = q^, defines our ariq^)- Numerically, 

ar(M^) = a.(M^) [l+0.5570a,(At^) + 1.702a,2(At^)] . (1.3) 

This exact (by definition) relation is used to run aj- between different scales using the intermediate 
a, coupling, whose running can be reliably performed at 4-loops over the range of scales relevant 
to the observables considered. The HPQCD coupling, ay, is defined as follows. Beginning with 
Eq. (pr2|), one takes the RHS, with jx^ = e^^/^q^, to define an intermediate coupling, ay{q^). This 
coupling has a /3 function, j8', with known values of jSg, ■ ■■ but also non-zero, but unknown, 
higher order coefficients, {5^^ , whose values depend on the presently unknown j84.5_.... The final 
HPQCD coupling, ay, is obtained from «(, by adding terms of 0(o;f ) and higher with coefficients 
chosen in such a way as to make j8j = jSj^ = • • • = 0. Since j84,5 ... are not known, the values of the 
coefficients needed to implement these constraints are also not known. 

Using the expansion parameter ar, no perturbative uncertainty is encountered in converting 
the fitted reference scale ar value to the equivalent reference scale MS result. This is not true 
for the HPQCD parameter ay. Higher order perturbative uncertainties, however, do remain in our 
analysis. To see where these occur, and to understand the motivation for the alternate HPQCD 
choice, let us define ao = &{Qq), where Qq is the maximum of the BLM scales (corresponding 
to the finest lattice) for the observable in question. We next expand the couplings at lower BLM 
scales (coarser lattices) for the same observable, in the standard manner as a power series in ao, 

cc{Qk) = £ PN{tk)a^ (1-4) 

where tk = log {qI/Qq), and the pN{t) are polynomials in t with coefficients determined by those 
of p. Substituting this representation into Eq. ( |1.1| ), one obtains the following expression, where 
we replace any occurences of j8o, ■ • • , j82 with their known numerical values and display only those 
terms involving one or more of and the unknown quantities j84,j85, • • ■, C3 , • • •: 

^ = • • • + (4 (^cf + •••)+ tto^ (^cf - 0.01027j83 - 2.865cf + •••)+ (^cf - 0.00321 p^tk 
-3.58lcfhk + [0.02573tl-0.02053cf\k]p3[5.l29tl - 1.621f;t]cf + • 
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+aj (^cf^ -0.00m0p5tk + [0.00936ltl -0.006536cfhk]p4 
+ [-0.04213?! + (0.01664 + 0.06617cf* )ff -0.03080cf ?/t]i83 -4.297cff*: 
+ [7.694ff - 2.026fi]cf + [-7.347f| + 6.386ff - 43S2tk]cf^ + •••) + ••■. (1.5) 

Running the MS coupling numerically using the 4-loop-tmncated jS function is equivalent to keep- 
ing terms involving jSo, • • • , jSs to all orders, and setting j84 = jSs = • ■ • = 0. The neglect of j84, jSj ■ • • 
also alters j34,j85,- -, and hence produces a "distortion" of the true f<;-dependence, beginning at 
0(0:0). Since it is the scale-dependence of Ok which allows one to fit the unknown coefficients 
cf\ as well as Oo, it follows that the 4-loop truncation forces compensating changes in at least 
the coefficients c^l^l . A shift in the values of cf\ however, leads also to a shift in the C?(o;q) 
coefficient which, in general, will necessitate a compensating shift in cf^ as well. This in turn will 
necessitate a shift in ao- Since the MS j8 function is known only to 4-loop order, such truncated- 
running effects are unavoidable at some level. From Eq. ([L^), however, it follows that their size can 
be minimized by taking Qq as large as possible (achieved by working with the observable with the 
highest intrinsic BLM scale) and keeping from becoming too large (achieved by restricting one's 
attention, if possible, to a subset of finer lattices). Note that, by defining ay in such a way that the 
4-loop-truncated p function is exact, the HPQCD expansion parameter choice, by definition, 
avoids these truncated-running problems. The price paid is the unknown relation between tty and 
a, beyond 0{aj). The impact of this uncertainty on as{Mz) cannot be controlled, either through 
the choice of observable or through the restriction to a subset of finer lattices. 

From the discussion above we see that the two different coupling choices lead to complemen- 
tary analyses. If the impact of the neglect of higher order perturbative corrections in both cases is 
small, the two approaches should give compatible results for analyses based on the same observ- 



ables, providing a form of mutual cross-check. Good agreement is, indeed, found (see Ref. [ |15| ] 
for details). 



2. Results 

We now turn to the results of our analysis. We employ data on our observables from the MILC 



a ~ 0.06, 0.09, 0.12, 0.15 and 0.18 fm ensembles []19|]. To minimize incompletely incorporated 
higher order perturbative contributions, we perform our main analysis using the three finest lattices, 
expanding to full 5-fold fits to test the stability of our solutions. The physical scales for the various 
ensembles are determined using the measured values of r\/a and the recent MILC assessment. 



ri = 0.318(7) [20]. The uncertainties on the extracted a, values associated with those on ri/a and 
ri are added linearly to arrive at an "overall scale uncertainty" contribution to the total error. 

Quark-mass-dependent non-perturbative contributions are found, using the data for different 
mass combinations ann /anis, to be very linear in lam^ + anis, allowing these contributions to be fit- 
ted and removed with good reliability. Because of details of the analysis not discussed here, the un- 
certainty in this subtraction appears as part of "the overall scale uncertainty" discussed above [jl^]. 

Mass-independent non-perturbative contributions are assumed dominated by the D = 4 gluon 
condensate contribution. The corresponding leading order contribution to the m x n Wilson loop 
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Wmn, denoted by SgW^n, is known from Ref. pT|], 



2 

~7l 2 2 4/ /-<2 \ 



8gW,„n = 36-"^ « « ) • (2.1) 

As central input for the condensate we employ the result of the updated charmonium sum rule 
analysis of Ref. {^G^) = (0.009 ± 0.007) GeV"^. Since the error is already close to 100%, we 



take the difference between results obtained with and without the resulting subtraction as a measure 
of the associated uncertainty. The neglect of mass-independent non-perturbative contributions with 
D > 4, for which no pre-existing constraints are available, should be safe so long as the estimated 
D = 4 gluon condensate subtraction is small. 

The shifts associated with the central gluon condensate subtraction grow with increasing lat- 
tice spacing a and, for the observables considered here, are (i) —0.01%, —0.02% and —0.08% 
for a ~ 0.06 fm, (ii) -0.1%, -0.4% and -1.3% for a ~ 0.12 fm, and (iii) -0.5%, -1.8% and 
—5.6% for a ~ 0.18 fm, where in each case the values quoted correspond to the observables 
log(Wii), log(Wi2) and log^Wn/ufi), in that order. The corrections, as claimed, are small, making 
the mass-independent non-perturbative subtraction safe for both the central 3-fold and extended 5- 



fold fits [23]. The subtractions are particularly small for the three finest lattices and for the plaquette 
observable, log(Wii). 

In line with the results of Ref. |^], we find that, even for the highest-scale observables and three 
finest lattices, the known terms in the perturbative expansion of the Ok are insufficient to provide 

(k) 

a description of the observed scale-dependence. When C3 is added to the fit, however, very good 
fit qualities are found, with x^/dof < 1 (very significantly so for the 3-fold fits). It follows that, 
with our current errors, it is not possible to sensibly fit additional coefficients in the expansions of 
the Ok- This raises concerns about possible associated truncation uncertainties. Since the relative 
weight of higher order relative to lower order terms grows with decreasing scale, the comparison 
of the results of the 3-fold and 5-fold fits provides one handle on such a truncation uncertainty. 
If higher order terms which have been neglected are in fact not negligible, then the growth with 
decreasing scale of the resulting fractional error should show up as an instability in the values of the 
parameters extracted using the different fits. We see no signs for such an instability within the errors 
of our fits, but nonetheless include a component equal to the difference of central values obtained 
from the 3-fold and 5-fold fits as part of our error estimate. This "stability component" is added in 

quadrature with the overall scale uncertainty, the gluon condensate subtraction uncertainty, and the 

(k) (k) 

small uncertainty associated with varying the C2 (and, if relevant, c\ ) within the errors in their 
numerical evaluations to arrive at the total error on our results. 

After converting our result for the reference scale = 3 coupling to the corresponding 
nf = 3 MS value, we run the result to Mz using the usual self-consistent combination of 4-loop 
running and 3-loop matching [^4|], taking the flavor thresholds to lie at rmc{mc) and rmb{mb), with 
r varying between 1 and 3, mdmc) = 1.286 ±0.013 GeV and mb{nih) = 4. 164 ±0.025 GeV 
The evolution to Mz produces an additional 0.0003 uncertainty on as{Mz) 

Our results for a^(Mz), combining all errors as indicated above, are 0.1192(11) from the fit 
based on log(Wii), and 0.1193(11) from those based on log(Wi2) and log(Wi2/Mo)- The consis- 
tency represents an improvement over that of the corresponding results quoted in Ref. [|3|]. In line 
with the arguments above, we believe the most rehable analysis to be that obtained using the three 
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Figure 1: Contributions to the errors on as{Mz)- Shown are the results for as{Mz) obtained using (i) the 
3-fold fit strategy, with central values for all input (the "central" case), (ii) the alternate 5-fold fit strategy, 
still with central values for all input, and (iii) the 3-fold fit strategy, with, one at a time, each of the input 
quantities shifted from its central values by IcJ, retaining central values for the remaining inputs. The error 
bars shown in each such case are those associated with the uncertainties in ri/a for the various ensembles. 



finest lattices, and the observable, log(Wii), having both the highest intrinsic scale and smallest 
gluon condensate subtractions. Our final determination, based on this case, is thus 



a,(Mz) =0.1192(11) 



(2.2) 



The dominant contribution, 0.0009, to the total error is that associated with the overall scale uncer- 
tainty. A graphical depiction of the various components of the error is given in Figure |l| where, for 
clarity, only one-sided errors are shown. The results are in excellent agreement with those of the 
independent determinations mentioned above, whose average is shown in the figure by the shaded 
band. 
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